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Abstract 

A version of the Swiss-cheese model is investigated. The flat Friedmann- 
Robertson- Walker (FRW) universe is modified by the addition of several 
spherical regions with Lemaitre-Tolman-Bondi metric. 

We discuss light propagation in this model in detail to pave the way 
for a detailed numerical study of the Hubble diagram. 

1 Introduction 

Type la supernova data interpreted in terms of homogeneous, isotropic cosmo- 
logical models necessitate the introduction of a huge amount of dark energy, 
i.e., a substance of negative pressure (like the cosmological constant). A re- 
cent alternative suggestion is the consideration of inhomogenity in the density 
distribution. Indeed, lensing effects of local matter abundances modify the posi- 
tions of distant objects on the Hubble diagram. There are several recent papers 
about this question without a concordance whether this effect is strong enough 
to give account of supernova data,[] albeit the majority of the authors seems 
to deny this possibility. The question is rather important since the existence 
of dark energy (or cosmological constant) profoundly modifies our knowledge 
about matter, already challenged by the existence of dark matter. 

In the present paper we consider the problem within the framework of the 
exactly solvable Swiss-cheese model. The version of the model we use is con- 
structed the following way. 

Nonoverlapping spheres are cut from a flat Priedmann-Robertson- Walker 
(FRW) universe. The mass they contained before is compressed within each 
sphere to a smaller sphere with homogeneous density distribution. Hence the 
inner spheres form sections of some closed FRW model. Between the outer 
and inncir spheres there is a vacinim, where, due to spherical symmetry, the 
Schwarzschield metric describes the gravitational field. Within the inner spheres 
the closed FRW metric is valid, while outside the cut spheres the flat FRW 
metric is relevant. The metric and its first derivatives are continuous across the 
bordering surfaces of the different regions. 



We use the Landau conventions, i.e., wc assume H signature for the 

metric. The zeroth component of a four-vector is timehke, the first, second and 
third components are spacelike. Four-vectors are indexed by Latin letters, three- 
vectors by Greek letters. Throughoiit wc use c = 1 units. At light propagation 
we use the index 1 for labeling initial quantities and the index for labeling 
final (i.e., present) quantities. 



2 Light propagation in flat FRW universe 

Light propagation in gravitational field is governed by the Sachs optical equa- 
tions reviewed below. 

A light ray may be parametrized along its path by a parameter A so that its 
four- velocity is defined by 

It is a null vector, 

u'ui = (2.2) 
and satisfies the geodesic equation 

= . (2.3) 

Let us introduce another independent null vector along the path. It satisfies 

w*Wi = , (2.4) 

w'ui = 1 (2.5) 

and 

w'^^u'' = . (2.6) 

Furthermore, let us define two spacelike unit vectors {oi,P = 1,2) that are 
orthogonal to both the null vectors and w': 

Lim = (2.7) 

i>i = (2.8) 

LlL.p = -S^p (2.9) 
We also require that the vectors are parallel translated along the light ray: 

LU^" = (2.10) 
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Eqs. ( plol ) ensure that if Eqs. ([2j]), (|2^ are 

satisfied at one single point of the Ught ray, they wiU be satisfied at any other 
point of it as well. For the separation f * of two light rays we assume 



and 



hence we have 



Cu, = 



Cm = , 



Q = l,2 



(2.11) 
(2.12) 
(2.13) 



The coefficients da describe the proper separation of the two nearby light rays. 
For their derivative with respect to the path parameter A we have in the most 
general case (since no rotation is possible) 



dda 
~dX 



0=1,2 



where Wa/3 is symmetric and traceless: 



P o- 

cr -p 



Now, the geodesic deviation equation 



(2.14) 



(2.15) 



(2.16) 



leads to the Sachs optical equations for the expansion rate Q and shear Wap'. 



dA+^ +2- 



(2.17) 



dX 

where 



2ewafi - C,jfc,L>%'=4 (2.18) 



W = WapWpa (2.19) 

and Cijki stands for the Weyl tensor. 

The physical situation (e.g., that the two nearby light rays are emitted from 
a point source) is specified by the initial conditions. After having solved these 
equations, Q governs the change of the cross section A of a light beam: 

d A 

^ = 2QA . (2.20) 
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In a flat FRW universe the metric can be written as 

= dt^ - a{tf {dr'^ + (di?^ + sin^ M^p'^)) . (2.21) 
In a matter dominated universe the scale factor a{t) may be given by 

2 

a{t) = a{to)(^lHoty (2.22) 

where 

H'o = ^p(io) (2.23) 

according to the first Friedmann equation, Hq denoting the present value of 
Hubble's parameter. 

The following formulae hold, provided that the light is propagating radially 
away from the origin. 

The four-velocity is given by (up to a constant multiplier, which determines 
the color of the light at emission) 

u* = - (2.24) 
a 

u"- = \ (2.25) 



,1? 



For the trajectory we get 



(2.26) 
«^ = (2.27) 



r = ^ (2.28) 



or 



r dt 

For the expansion rate we have 



a{t) a{ti) 



(2.29) 



e=4 + ^ = ^ (2.30) 
and thus the cross section area of the beam is 

A = OaV^ (2.31) 

The redshift is defined by 

1 = 4^_1 , (2.32) 



w*(to) a{ti) 
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For the intensity we have 



I = L 



where 



47rA(to)(l + 2)2 ' 



(2.33) 



(2.34) 



and L stands for the absolute luminosity. Then the luminosity distance is given 

by 



dL 



-\/rT^(x/rT^-i) . 



47rJ Hq 

This is the Hubble- diagram in a matter-dominated flat FRW universe. 



(2.35) 



3 4 5 

redshift 



Figure 1: The Hubble-diagram in a matter-dominated flat FRW universe. 



Let us describe light propagation in an arbitrary direction. As space is 
homogeneous and isotropic, it is enough to perform a translation of the origin. 
This means that 

r = r' + rfe (2.36) 
Suppose that in the coordinate system K' the light ray is radial. 
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The above coordinate transformation implies that 

r' = Y^r^ + r| — 2rrfc (cost? cos i?fe + sin i? sin i?^ cos{(p — (fk)) (2.37) 

n/ rcosi?-rfcCost?fe 
cosw = . (2.38) 
_^j.2_ 2rrk (cos i? cos i?fe + sin ■!? sin i?/; cos(i^ — Lpk)) 

, r sint^sina? — Tfe sini^fe sinccfc „„x 

tant^' = — (2.39) 

r sm y cos (p — rk sm cos (pfc 

The trajectory is given by the inverse coordinate transformation 

r = + ^fc + ^f'^k (cos I?' cos t?fc + sin -d' sin i?^ cos(</5' — ipk)) (2.40) 

r' cos + Tfe cos ??fe 

cos 17 = (2-41) 

■\/ + J^fc + 2r'r/c (cos z?' cos 'dk + sin i?' sin -dk cos(</3' — ifk)) 

r' sin 1?' sin + sin 1?^ sin y>fc 

tano? = -— — ; : — (2.42) 

r sm cos ip' + rk sin Vk cos ifk 



where i?' and if' are constant, while 



t 



a{t) a{ti) 



(2.43) 



The expansion rate is 



e=4 + 47 = ¥^ (2-44) 
and for the cross section we have 

A = na^r''^ (2.45) 

Let us determine the transformation of the velocity. The timelike component 
is clraly unchanged. Covariant spatial vectorcomponents transform as 

vo. = -^v'p (2.46) 
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The transformation matrix elements are given explicitly by 



dr 

dip 
dr 



dip 



dip' 

dr 
dip' 
'&& 

dip' 

dip 



r — rk (cos "9 cos + sin sin -(5^ cos ((/? — ipk)) 
-\/r^ + — 2rrfc (cos i9 cos i9fe + sin d sin dk cos(i^ — ^k)) 
rrk (sint?cosi?fc — cos sin t^^; cos(<^ — ipk)) 
_|_ ^2 _ 2rrk (cos 1? cos iJfe + sin i? sin 1?^ cos{ip — tpk)) 
rrk sint^sini?/; sm{ip — ipk) 
^^2 _|_ ^2 _ 2^^^ (qos COS + sin ?9 sin i?^ cos((p — </?fc)) 

— (r cos'i?fc + Tfc cos t?) + (r cos i9 + rk cos -dk) (cos?9 cosi^j. + sin?9 simJ^ cos((^ — <^fc)) 
(r^ + — 2rrfe (cost? cos + sinz?sint9;j cos{ip — ifk))) 



■ sin^ ^ + sin^ — 2rrfe sin i? sin i?/; cos((^ — yi/s) 



sin I? + 



r/s (r C0S1? — rfc cos??fe) (sint?cosi?fe — cost? sin ■iJfc cos(<p — ipk)) 

(r^ + — 2rrfe (cos i? cos + sin i? sin i?/; cos(y — ipk))) 
r 

X y 

\/ sin^ ?9 + sin^ i?^ — 2rrfe sin i9 sin t^fe cos((^ ^ '■Pk) 

r (r cost? — rk cosdk) sint^sinz^fe sin((p — ipk) 
(r^ + — 2rrfc (cost? cost?/; + sini?sint?fe cos(<^ — i^fe))) 

J sin^ ^ + ^fc sin^ t?fe — 2rrfe sin t? sin dk cos((^ — ipk) 
rk sin 1? sin t?fe sin((^ ~ Vk) 
sin^ + '''fc sin^ i?fc — 2rrfe sin t? sin t?^ cos((p — ipk) 

rrk cos t? sin t?fc sin((^ — (^fe) 
sin^ 1} + rl. sin^ t?^ — 2rrfe sin t? sin t?fe cos((^ — ipk) 
sin^ 1? — rrk sin t? sin t?fe cos((^ — (^fe) 
sin^ + r1 sin^ t?fe — 2rrk sin i? sin t?fe cos(<^ — Vk) 

The transformation of contravariant vectors may be derived by pulling in- 
dices down, transforming the covariant vector thus obtained, and finally pulling 
indices up. Note that the metric is the same both before and after the transfor- 
mation, if expressed in terms of the corresponding coordinates. Hence for the 
four- velocity we get 

1 



u = 



uf = 



r — rk (cost? cos t?fe + sin i? sin t?^ cos((y9 — ipk)) 
g2^y.2 2rrk (cos i? cos 1?^: + sin t? sint^fe cos{ip — ipk)) 

rk (sin t? cos t?fc — cos t? sin t?^ cos((/3 — ipk)) 
a^r^Jr"^ + r| — 2rrfc (cost? cos i?fe + sin t? sin i?^ cos(</3 — ipk)) 

V]; sill 0]. siii(^ — y /.) 
a^r sin ■d^Jr'^ + ''fc ~ ^rrfc (cost? cos t?fc + sin i? sin t?^ cos((^ — ipk)) 



(2.47) 
(2.48) 

(2.49) 

(2.50) 
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3 Light propagation in closed FRW universe 

In a closed FRW universe the metric can be written as 

ds^ = de - a{t f (^^^ + {dd^ + sin2 M^^)^ . (3.1) 

In a matter dominated universe the scale factor a{t) may be given parametri- 
cally: 

a{t) = ^sin^x (3.2) 
t = /5(x-sinxcosx) (3.3) 



where 



P=^-^p{to)a\u) . (3.4) 



Note that according to the first Friedmann equation 

P ..2,1 



q/, \ — ^ „, , , (3.5) 

Hq denoting the present value of Hubble's parameter (which may deviate 
from that in the flat FRW region). Further, we have 



X=17W^ ■ (3.6) 
2/9 sm X 

Provided that the light propagates radially in the chosen coordinate system, 
d'd = Q and dip = holds, hence ds^ = implies 

dt^ = a{tf . (3.7) 

The solution is 

r = sin(2(x-Xi)) , (3-8) 

whore xi = x(ii)- 

The components of the four velocity are (up to a constant multiplier again) 

w* = ^ (3.9) 

vT = ^ / (3.10) 

= Q (3.11) 

u'P = Q (3.12) 
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The cross section area of the hght beam emanated from the origin may be 
expressed as 



which imphes according to Eq.(2.20l that 



e 



(3.13) 



(3.14) 



Then 



e + ae^ = 



3 /3 



i.e., it satisfies Sach's optical equations (2.17), (2.181. 
The redshift is given by 



^ _ a(^o) _ _ sin (xo) 



a(ii) 

For the himinosity distance we have 



sin^(Xi) 



1 



(3.15) 



(3.16) 



rfL = (l + ^)\/^=(l + ^)a(to)ro = /3^!4M^i"(2(xo-Xi)) ■ (3.17) 
V " sm (xi) 

Eqs.( |3.16| , (3.17) express the Hubble diagram in parametric form, the parameter 
being xi- Exphcitly, using Eqs.(3.2|, (3.5), (3.8), too, we get for the expanding 
case (xo < 7''/2) 



gga§((l + z)(l + gg«^) 



Hpao (1 - H^al) + (1 + z)H^al 



(3.18) 
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redshift 

Figure 2: The Hubble diagram in a closed, expanding FRW universe. Different 
curves correspond to different ifooo values. Note that in the limit Hoao 1 
the Hubble diagram of the flat FRW is obtained. 

In the collapsing case (tt > xo > 7i'/2) the parametric form is more suitable, 
since the Hubble diagram becomes multivalued, because light may reach the 
observer from the opposite direction, through the rest of the universe as well. 
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redshiit 



Figure 3: The Hubble diagram in a closed, collapsing FRW universe. Different 
curves correspond to different HqUo values. Note that the redshift can be nega- 
tive (i.e., we may get a blucshift), and in that region the diagram is multivalued. 

We also need the expressions for light propagation through an arbitrary point 
in an arbitrary direction. Due to homogenity and isotropy, this can be obtained 
from the expressions above by performing a suitable rotation of the coordinate 
system (which is an isometry). This transformation may be represented as 
rotation of a three-sphere embedded into a four dimensional Eucledian space. 
Let the Cartesian coordinates of this space denoted by xi, X2, X3 and x^. The 
homogeneous and isotropic three dimensional space of our interest is the surface 
of a sphere of radius a. This surface may be parametrized as 

xi = a sin sin cos (3.19) 

X2 = a sin ?/> sin sin (/? (3.20) 

X3 = a sin?/) cos?? (3-21) 

X4 = acos'0 (3.22) 

Obviously, the usual radial coordinate r is given by 

r — sin . 



The line element on the surface is just the spatial metric we get from Eq. (3.1) 

i.e. 

ds^ ^ + r\dH + sin2 M^^) 

\ 1 — 



(dV + sin^ '4}{d^d -t- sin^ M^Lp)) (3.23) 
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A four-dimensional Euclcdian rotation around any axis through the origin 
is an isometry of the three-surface. We shall use three such rotations in order 
to transfer the origin of the three-space (i.e., r = 0) to another point whose 
coordinates are V = V'fe) = "^k and ip = (fiki 

1. 2:3 = const. , Xi = const. 

This readily implies that r (or tj)) and 1? are constants, too, while changes 
according to 



cos \ _ I cos 77 — sm 77 \ / cos (/? 
sin / \ sin rj cos 77 j \ sin ip 



(3.24) 



where 77 stands for the angle of rotation. This implies 

= <^ + 77 , 

hence choosing 

77 = -<pk (3.25) 

sets the coordinate of the new (primed) origin to yfe. 

2. X2 = const. , Xi = const. 

This readily implies that r (or ij)) is constant, too, while 1? and </? change 



according to 

sin 7?' cos \ _ ( cos 7 — sin 7 \ / sin 7? cos ip 
cos 7?' J Y sin7 cos 7 J \ cost? 

where 7 stands for the angle of rotation. This implies 



(3.26) 



cos 7?' = sin 7 sin 7? COS 1^ -|- COS 7 COS 7? (3.27) 

, COS 7 sin 7? cos (jj — sin 7 cos 7? ,„ „„n 

cos if' = . ' y / (3 28) 



1 — (sin 7 sin 7? cos (p + cos 7 cos 7?) 
If (p is already zero (due to the previous transformation), this implies 

7?' = 7? - 7 , <p' = 0, 

hence choosing 

7 = 7?fe (3.29) 

sets the t? coordinate of the new (primed) origin to 7?^ while keeping ip 
unchanged. 
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3. xi = const. , X2 = const. 

This readily implies that ip is constant, too, while a? and ip change accord- 
ing to 



sin tp' cos ^' \ _ f cos v — sin \ / sin cos i? 
cosip' J y smu cosu J y cos'^ 

where v stands for the angle of rotation. This implies 



(3.30) 



cosV'' = sin sin COS + COS cos V' (3.31) 

cos i> sin -0 cos i? — sin ly cos ih 
cosy = — ^^^^^^^^^^^^^^^^^^^ (3.32) 



1 — (sin V sin ip cos ^ + cos cos ip) 
If 1? is already zero (due to the previous transformation), this implies 

Tp' = tp -u , ^' = , 

hence choosing 

V = ipk (3.33) 

sets the ip coordinate of the new (primed) origin to '(pk while keeping 
and ip unchanged. 

The succession of the above three rotations in the given order and with the 
given parameters corresponds to a spatial translation in the closed three-space. 
The resulting transformation may be expressed implicitly by 

sin ■0' sin cos = sini/» (— sini?fc cosi? -|- cosi^fc sin'!?cos(<^ — ly^fe)) (3.34) 
sin Ip' sin i?' sin ip' = sin ip sin -& sin(<y? — ^Pk) (3.35) 
sin V'' cos 1?' = sinV'cos^/^fe (cosi?/j cos^? + siniJfe sint?cos(<^ — i^fc)) 

— cos ^ sin -0^ (3.36) 
cosip' = sini/)fe sinf/; (cosi?fe cost? + sint?;; sini?cos(i^ — 99fe)) 

+ cos t/ifc cos (3.37) 

Similarly, the inverse transformation reads 

sin sin 1? cos y> = sini^fe cos (cos V-'fe sin-0' cos + sinV'fc cosi/'') 

+ sin tp' sm-d' {— sin ipk sin ip' + cos dk cos ipk cos </?') (3.38) 
sin ■)/) sin -i? sin = sint?^; sin (cos'f/'fe sin-)/)' cosi?' + sint/i^; cos t/;') 

+ sin tp' sin i?' (cos ipk sin ip' + cos i?fe sin c/j^ cos p') (3.39) 
sin tp cos 1? = cos 1?^ (sin ipk cos 'tp' + cos ipk sin ip' cos •&') 

— sm.ip'shi'd'cosif'smdk (3.40) 

cos Ip = cos Ipk cos tp' — sin ipk sin ^' cos i?' (3-41) 
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These expressions describe the trajectory of the hght when tf' and are 
constants and 



^' = 2{x-Xi) 
t = /3(x ~ sinxcosx) 

i.e., when the hght ray moves radially in the K' frame. 
The expansion rate is 



(3.42) 
(3.43) 



e = 



a 



\J\ — r'2 a + cot V'' 



a? r' 

and for the cross section we have 



A = na^ r'^ = Vla^ sin^ 



(3.44) 



(3.45) 



Finally, we calculate the components of the four-velocity of the light. We 
use the same technique, as in the case of the flat FRW universe. In the K' frame 
we have 



1 

a 

cos "ip' 

-.2 



a" 



,i>' 



Pulling down indices: 



Uf 
Ur' 



COS tp' 



Applying Eq.(2.46) we get in the K frame 
1 

ut = - 



(3.46) 

(3.47) 

(3.48) 
(3.49) 

(3.50) 

(3.51) 

(3.52) 
(3.53) 

(3.54) 



cos-iAfc tani/j — sinipk (cosi?/t cost) + siv.'&k smi!)coa(ip — ipk)) ,„ rr\ 
Ur = : — (3.55) 



u,„ = - 



sin ip' 

sin ipk sin tp (cos sin — sin i?^ cos "i? cos((^ — tpk)) 
sintp' 

sin V'fe sin tp sin sin sin(iy9 — ipi^) 
sintp' 



(3.56) 
(3.57) 
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Finally, by pulling up the indices with the contravariant metric tensor we get 
^ (3.58) 
cos V'fc sin xj) 



u = 



u = 



cos V' 



a? sin ij:' 

— sin cos V' (cos t?fe cos i? + sin 'd^ sin •& cos(<^ — Vk)) (3.59) 

sin il^k (cos iJfc sin i? — sin 1?^ cos 'd cos((^ — ^k)) 

sin i/j' sin ■)/' 
sin sin t?/c sin((^ - ipk) 



a? sin -0' sin tj) sin i? 



(3.60) 
(3.61) 



4 Light propagation in the Schwarzschild metric 



The Schwarzschild metric can be written as 



r / 1-'- 
Hence the metric tensor is given by 



sm' 



9u = 1-^ 
r 



9rr — 



1 



The non-vanishing components of the Christoffel symbol read 



r* 

tr 



-pr _ ' g 



pi? 

1 

■pip 



rg{r-rg) 

rg-r 

{Vg — r) sin^ 

— sin 1? cos I? 



vd (pr ri^ 



coti? 



(4.1) 

(4.2) 
(4.3) 

(4.4) 
(4.5) 

(4.6) 

(4.7) 

(4.8) 
(4.9) 
(4.10) 

(4.11) 
(4.12) 



Since we are dealing here with a vacuum solution of Einstein's equations, Ricci's 
tensor vanishes, hence Riemann's tensor is equal to Weyl's tensor. The nonzero 
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components are 



Rtrtr — Ctrtr — ^ (4-13) 

rg{r-rg) 

2r2 



Rm'& = Cfdfd = ~ ^ o 2 ^ (4-14) 



■n-tiptip — '^t<pt(p — {^■^'^) 



RrM — — ^,7 ^ (4.16) 

2(r - rg) 



r„ sin^ ■& 

7) 

R'&<P'&<fi = C^^^^ = -TgrsiT? -& (4-18) 



Rriprip — Criprif ~ 2{r — T ) ('^•l'^) 



where 



Rijkl — —Rjikl — —Rijlk (4.19) 

Rijkl = Rklij (4.20) 

In a spherically symmetric field particle trajectories (including photons) are 
confined to a plane perpendicular to the conserved angular momentum vector. 
Provided that the angular momentum vector points towards the z-axis and the 
center of the field is at the origin, = f during the whole motion. The motion is 
fully integrable. Integration may be performed by solving the Hamilton-Jacobi 
equation which describes light propagation if m = 0. Let us introduce the 
impact parameter p = ^ instead of the angular momentum J. The orbit is 
defined by 



which will be converted to an elliptic integral. To this end let us introduce 

(4.22) 



r 



^ = -1^^^=^ ■ (4-23) 



In terms of this new variable of integration we have 

The third order polynomial under the square root factorizes as 

?'-e' + 4 = (e-6)(?-6)(?-6) , (4.24) 
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where , ^2 and f 3 are the roots of the polynomial, namely, 



6 

6 



1 / ^ (a + 2iT 



1 + 2 COS 
1 + 2 cos 



a ~ 2tt 



where 



Now we can see that 



3^3 1 



2 p 



6 < < 6 < 6 < 1 



(4.25) 
(4.26) 
(4.27) 

(4.28) 
(4.29) 



Furthermore, during the motion < ^ < ^2 holds and the turning point of the 
light trajectory (i.e., its nearest point to the center of the fi eld) is at ^ = ^2- 
Note that for p < 3-\/3rg/2 no turning point exists (cf. Eq.(4.28)), thus light 
beams with such a small impact parameter inevitably fall into the center of the 
field. 

If the angle is measured from the turning point, 
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V(e-ei)(e-6)(e-e3) 



(4.30) 



Performing the integration we have 



(p = 



K 



™(t) 



-F 



'V3 e-Ci 



2 sin (^1) ' y sin I 



3 ) 



(4.31) 



Especially, when the light comes from the infinity, the change of the azimuthal 
angle until reaching the turning point is 



V2^ 



sm(^^) 



K 



sm I 



-F 




\/3 a 
1-^tan-, 



sin (5i±^) 



(4.32) 
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In the above equations F{z, k) stands for the incomplete elUptic integral of the 
first kind and K{k) for the complete elliptic integral of the first kind, i.e. 



K{k) 
F{z,k) 



F(l,fc) 



dt 



V(l-t2)(l-fc2t2) 



(4.33) 
(4.34) 



Gravitational field bends the light beam. Let us denote by Sip the deviation 
of the outgoing (to infinity) light beam from the original (ingoing) direction. 
Explicitly, we have 



6ip = 2ip — TT 



(4.35) 



The series expansion of Eq.(4.32 ) in terms of Tg/p yields the classical expression 
for the deviation of light. The zeroth order vanishes, 

<5<^(o) = 2^(°) - TT = , 

and in first order we get 



(4.36) 



P 

The time dependence r — r{t) is determined by 

dr 

ct 



(4.37) 



(4.38) 



or, using the new variable of integration introduced in Eq.(4.22)) 



ct p 



d^ 



(4.39) 



where the roots ^i, ^2 and ^3 are determined by Eqs.( 4.25 ^4.27 1. Performing 
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the integration we get 




(4.40) 



Here the + sign holds if the Ught is going away and the — sign, if it is approach- 
ing. The symbol C is a constant of integration. We choose it so that t — 
corresponds to the turning point, i.e., ^ = ^2- Hence 



C 



If /Q! + 27r 



2 . (a + ^\ 



K 



sin(f) 



sm 



6 



1 + 2 cos 



-H 



1 — cos 



' Q+27r 



-H 



-2y3sin(|) 
^l + 2cos(^i^^' 

V3sin(f) [ 

1 — cos 



sm 



' Q + TT 



(4.41) 



In Eq.(4.40 ) E(^z^ k) stands for the incomplete elliptic integral of the second kind 



and n(z, k) for the incomplete elliptic integral of the third kind. I.e. 
E{z,k) = / \l \ dt 



l-t2 



z 

Ii{z,v,k) = J 



dt 



(1 - I^t^) y^{l - t^){l - kH^) 



(4.42) 
(4.43) 
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In Eq.(4.41) E{k) stands for the complete elliptic integral of the second kind 
and n(i^, k) for the complete elliptic integral of the third kind, i.e., 



Eik) = E{i,k)= I 



1 - /c2i2 



-dt 



\l{v, k) 



i 

n(i,;.,fc) = j - 



dt 



(1 - Vt^) ^(l-<2)(l_fc2t2) 



Eqs. (4.21), (4.38) imply 





dr 




dip 




dr 




~ ~dt 


7r 

~ 2 ' 


we have 



y^p2 ( r ) 



= 1- 



1 - 



1-^ 

r 



u''=0. 

Therefore, the r component of the geodetical equation sounds 



dvJ 
dr 



which implies according to Eqs.(|4.6 - 4.9) 

1 r, 



du" 
dr 





u"' 2r2 ^ r ' 



1 



2r^ 



(1 



Using Eqs.( |4.46[|4.47[ ) we get 

dr 

with the solution 



Mo CXp 



dr 



Wo \/l - ^ ( 1 - 



(4.44) 
(4.45) 



(4.46) 
(4.47) 

(4.48) 

(4.49) 

(4.50) 
(4.51) 



(4.52) 



where Uq is an integration constant. Consequently, the four-velocity of the light 
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in Schwarzschild metric is given by 



M = uo ^ (4.53) 



= no ^1-^(1-^) (4.54) 
= (4.55) 
= Mo ^ (4.56) 

In order to determine the light intensity during propagation one has to calculate 
the change of the cross section of the light beam. Suppose that the centerline 
of the narrow beam of an elliptic cross section lies in the i9 = ^ plane. This 
plane is clearly a symmetry plane of the beam, hence one of the axes of the cross 
section ellipse remains collateral with this plane and the other is perpendicular 
to it (assuming initial conditions already having possessed this property). We 
have to consider two type of trajectories on the outer surface of the beam: 

• trajectories lying in the = f plane (type 1) 

• trajectories at a maximal distance from the = f plane (type 2) 
Certainly we have no further linearly independent possibilities for the deviations. 
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type 2 



symmetry plane 



type 1 



Figure 5: Cross section of the beam. The symmetry plane and orbits type 1 
and 2 are shown. 



Trajectories of type 1 lie in the same plane as the centerline, but their initial 
conditions and impact parameter slightly deviate from those of the reference 
centerline trajectory. The perpendicular distance dl from the centerline may be 
expressed as 



dl = r\ 1 



1 - 



(4.57) 



Similarly, at the initial point (we consider ip to be the function of r, p and the 
initial position ri and ipi) we have 



dli=riJl 



1 



dip dtp 
- — dpi + -TT—dri 
oipi ori 



(4.58) 



Also, we need to take into account that the total change in the value of the 
angle may be written as 



dp ~ -^dr 
or 



dip dip dip 
-T^dp + - — dipi + ^— ari 
op otp>\ or I 



(4.59) 



Combining Eqs. (4.571, (4.58), (4.59), (4.21) we get for the expansion rate for 
trajectories in the plane = f 



d[^ 



r\ I 



p- 



dp 
dh 



dr 



1- 



(4.60) 
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As for trajectories type 2, they also lie in some plane which makes some small 
angle d'y with the original = f plane. Simple geometrical considerations imply 
that the perpendicular distance from the reference trajectory may be written as 



r |sin((^ — $)| ^7 



(4.61) 



where $ denotes the azimuthal angle of the intersection line of the planes. Hence 
the expansion rate for trajectories not in the plane ?9 = ^ we get 



r |sin((^ — $)| 
ri |sin(^i - $)| 



(4.62) 



Note that eventually we let the initial point coincide with the light source that 
implies ^ = ipi. A direct substitution at this stage is obviously not yet allowed. 

Putting things together, if the light beam has cross section Ai at the point 
(ri , iy9i ) , then at another point (r, tp) its cross section is 



|sin(v3 - $)| 
ri |sin(v3i - $)| 



dp 
dh 



dr 



(4.63) 



Below we cast this equation to a more suitable form. Let us denote by fii 
the angle of the light ray and the radial direction at the initial point (ri, ipi) in 
the locally Eucledian frame. 

Geometry implies 

dri 



cos ^1 oc y/—grrdr = 



sin/ii c>c ridipi , 



1 - 



(4.64) 



(4.65) 



the proportionality constant being the same in both equations. Using Eq.( 4.2l"] ) 
we obtain 



cot ^1 
therefore 



dri 



dip I 



1 _ i V ^ 

ri 



1 - 



ri 



sm fXi 



P_ 
ri 



1 



LI 
ri 



(4.66) 



(4.67) 



Provided that the initial point is at a small distance L from the source we may 
express the aperture d/ii ^ 1 of the light as 

dh 



dpi = 



L 



(4.68) 
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Eqs.(|4.67|) and (|4.68|) readily imply that 

ri cos /ii 



dp 



Ljl 



(4.69) 



Suppose for simplicity that the light beam has a circular cross section near its 
source. Then the distance (4.61) may also be expressed as 



ri |sin((y9i — ^)\d'y = L sin /^id7 



(4.70) 



Near the source, when i — 0, in Eq.(4.63) the term containing dp/dli becomes 
dominant, so the other term may be neglected. Introducing the solid angle 



L2 



(4.71) 



Eqs.(4.63l, (4.69) and (4.701 lead to 



A 



r\rl |sin(^o - ^i)| ^^^1-^(1-^)^1-^(1-^) 



P 1- 



dr 



(4.72) 



Note that at this stage we have been allowed to let the initial point coincide 
with the source, hence ^ — (p\ has been substituted. Also, the final point is 
denoted by (ro, (/3o) henceforth. Using the equation of the path (Eq.( |4.31 )) we 
have 



Vsin((a + ^)/3) 



F 



sin {a/3) 



F 



The redshift is given by 



2 sin(a/3) ' V sin ((a + 7r)/3) 



^gtt{ti)u\ti) 



2 sin (a/3) ' V sin ((a + 7r)/3) 

(4.73) 



sin (a/3) 



(4.74) 



hence inserting u from Eq.(4.53) we get 



1 - 



(4.75) 
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Solving this for ri we have 



ri = 



1 



(4.76) 



The luminosity distance is 



(4.77) 



which, combining with Eqs.(4.72), (4.73) and (4.76) yields the Hubble diagram. 




0.5 1 1.5 2 

redshiit 

Figure 6: The Hubble-diagram in Schwarzschild-metric. The luminosity dis- 
tance is given in Vg units. 



In case of a different p we get a different Hubble-diagram. In Fig. ([4]) the lower 
envelope is determined by p = 0, which corresponds to radial light propagation. 
In that case we have 



which implies according to Eq.( 4.72| ) 



(4.78) 



(4.79) 
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which, together with Eqs.(4.76), (4.77) yields the luminosity distance 



ro 



(4.80) 



Concluding this section, we calculate the expansion rate Q by using Eqs.(2.20) 

dA 



( |4.54[ ), ( |4.56[ ) and ( |4.72| ). We get 

dA ^ 



2A dX 2A [ dr 



P 



uo \ ^cot((^- (^i) + 



2 1 



P_ 

r 



dip 

^ 2r 2r) 

1_ p!fl_i) 



(4.81) 



5 Matching conditions 

Matching the different regions at the bordering spherical surfaces can be done 
explicitly in terms of suitable coordinate transformations. This transformation 
is specified - to the necessary degree - by the requirement of continuity of the 
metric and its first derivatives when crossing the bordering surface. Below 
we list the expressions obtained. For completeness, derivation is given in the 
Appendices. 



5.1 Matching the flat FRW metric to the Schwarzschield 
metric 

Metric outside of the sphere (fiat FRW metric) : 

ds^ = dt)-a {tff {dr) + r) [dd^ + sin^ Mip^)) (5.1) 
Metric inside of the sphere (Schwarzschield metric): 

ds^ = (l - ^) dtl ~ - rl (d^2 + sin^ (5.2) 

The bordering spherical surface lies at — Ri. 
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Coordinate transformation: 



9 

4 2 



Hrf - R^)^P^^ + {r, - R,r ^ (5.3) 



«or/ + (r/-i?i)^^ (5.4) 

= (5.5) 
fs = (5.6) 

The matching condition is the following: 

alRlH^ = rg (5.7) 

5.2 Matching closed FRW metric to the Schwarzschield 
metric 

Metric outside of the sphere (Schwarzschield metric): 

ds^ = (l - [!^) dtl - ^pri - rl {M^ + sin^ M^^) (5.8) 
Metric inside of the sphere (closed FRW metric): 

ds^ ^dtl~a {t,f [0-2 + rl {d^^ + sin^ ^V)) (5.9) 

The bordering spherical surface lies at rc — R2- 

In this case a(tc) fulfills Friedmann's equation whose parametric solution is 
given in Eqs.([3^-(|3!6l). 

It is convenient to express the transformation laws in terms of this parametriza- 
tion. Introducing the constant angle C by 

i?2==sinC (5.10) 
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we get 



(3 



cos ^(1 + 2 sin^ Ox — cos C sin x cos x — 2 sin'^ C artanh 



, . , tan C sin"^ x cos x 

+ (rc - sm C ^- ^ — 

cos^ C ~ cos^ X 

1 . sin^ X cos X (cos'* C - cos^ x) 
+ -(?'c — smQ ^ 

2 cos'^ C (cos^ C ~ cos^ x) 

■ 2 3 , . 2 sinC 
r.sm x+^(r.-smC) ^ 



tan X 
^ tan C, 



(5.11) 

(5.12) 

(5.13) 
(5.14) 



Further, the consistency condition 



(5.15) 



must be fulfilled. The constant (3 is defined by Eq.(3.4). 



6 Light propagation and beam expansion in the 
Swiss-cheese model 

Now we put everything together: using the previous analytic expressions we 
describe how to determine redshift and light intensity if the light ray passes 
through several regions, each having one of the three metrics considered. Above 
we considered these regions separately, therefore, we have expressions for light 
trajectories, cross sections and expansion rates when the light source is in the 
same region. 



The next problem is what happens at the borders. There Eqs.(5.3 1-(5.6 1 or 



Eqs.(5.11 1-(5.14|, respectively, describe the coordinate transformation between 
the regions preserving the continuity of the metric and its first derivatives at 
the border. The four velocity transforms accordingly, namely. 



1 - 

ao 



U% 3 



43l fk.\ ' , 

to \to J 



(6.1) 



'^aoTftJtj^'u'f + aoi^-^^'u} (6.2) 

it = uj (6.3) 
it = u^f (6.4) 
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at the flat FRW-Schwarzschild border and 



uf 



t COS C sin^ X tanCsin^xcosx 

. 2 • 2 A ^cP .2 ■ 2 /- 

sm X ~ C x ~ sm 

u*rc cot X + sin^ x 



(6.5) 

(6.6) 
(6.7) 
(6.8) 



at the closed FRW-Schwarzschild border. The expansion rate and the cross sec- 
tion are invariant at the border since they are directly measurable quantities. 
Position, velocity, expansion rate and cross section render possible the calcula- 
tion of these same quantities beyond the border, inside the other region at a 
later time. Below we describe how this can be achieved. 

If the light ray enters the flat FRW region at time t = tin then Eq.(2.44| 



determines r' at this instant of time, given the expansion rate O. The quantity 
r' obtained would be the coordinate distance from the source had the light ray 
travelled all the time in flat FRW metric, like in a flat FRW universe. Next, 
using Eq.(2.43l with t = tin determines ti, which would be the emission time 



from the source in a flat FRW universe. Given the cross section A at the entrance 



time, Eq.(2.45l determines the solid angle Q of the light ray at the imaginary 



source. The light trajectory is given by Eqs. ( 2.40 )-( 2.42) where the parameters 
rk, iSk, fk, and t/s' are determined from the position and velocity aX t — tin, 
using Eqs. (|2.40|)-(|2.42|) and Eqs. (|2.48|)-(|2.50|) together, where r' is the value 



determined from the expansion rate as discussed above. Note that while there 



seems to be more equations than variables, in the reality Eqs. (2.48l-(2.50l are 
not independent, since they satisfy 

((w^)' + {u^'f + sin^ ^ {u^^f^ - {u'f . 

If the light ray enters a closed FRW region the procedure is analogous to that 



used in the flat case. Given the expansion rate 8, at entrance time Eq.(3.44| 



determines r', or, equivalently, tp' , the coordinate distance of the imaginary 
source (which would be there if light had travelled all the time in a closed FRW 
universe). Next, Eq.(3.43l determines the emission time ii, while Eq.(3.45l 



yields the solid angle il of the light ray at the imaginary source. The light 
trajectory is given by Eqs. (3.38l-(3.43) where the parameters i/ifc, §k, 'fik, ^' 



and if' are determined from the position and velocity at entrance time, using 



Eqs. (3.38)-(3.41) and Eqs. (3.59)-(3.61) together. 



Finally, when the light ray enters a Schwarzschild region, Eqs. (4. 81 1 and 



(4.731 determine ri and tpi (for the given the expansion rate Q), i.e., the po- 
sition of the imaginary light source. The emission time ti may be found from 



Eq.(4.40). Eq.(4.72) determines the solid angle fl of the light ray near the 



source. Eqs. (4. 31), (4.40) describe the path and the time dependence, while 



Eqs. (4. 81 1, (4.72) yield the expansion rate and the cross section at any later 
time. 
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If there are several spheres ("holes in the Swiss-cheese"), the above proce- 
dures should be applied in succession. Since the matching formulae assume that 
the origins of the different coordinate systems coincide, shifting from one sphere 
to the other necessitates a spatial translation in the flat FRW region, given by 
Eqs.( |2.37| )-( |2.39| ) or ( |2.40[ )-( |2.42| . Here r^, and (pfc are the coordinates of 
the centre of the next sphere in the old (unprimed) coordinate system. The 
velocity components transform accordingly. Further, a rotation around an axis 
through the origin may also be necessary afterwards, in order to ensure that 
the light trajectory lies in the plane = 7r/2, this being the situation when our 
expressions for the Schwarzschild region hold true. 

Finally, we get at the observation point A{to) for the cross section and 
Ut{to) / \/goo(to) for the frequency. The redshift is given then by 



z = , — i = , — i 



and the luminosity distance by 



= (1 + z) 



Aito) 



n 



(6.9) 



(6.10) 



These expressions complete the determination of the Hubble diagram. 
A detailed numerical investigation based on the equations presented is under 
way. 
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A Matching the flat FRW metric to the Schwarzschild 
metric 



At the bordering spherical surface the line element is 









\dtf 


1 




1_ i 


\dtf 



, dts dts 
^dtf drf 



dti + 2^^dtfdrf 



dr ' I ^ 



<f 



-rf (rfi?2 „ sin^ ?9 (^.1) 
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and the metric is given by 



.910 



9oo 



901 



911 



922 = -r 



1-^ 



dtf 



drs 

1 - i V at 



= 1 



1 



1 - 



dts^ 
drf 



<9ts dts 

2 1 / n \ 2' 



rt=Ri 



1 



1 



(A.2) 
= (A.3) 
(A.4) 



933 
902 



= -r: sin -d 



s\rf=Ri 
2 „• 2 , 



tjRi 



rj=Rl 



520 — 903 — 930 — 



Eqs.(A.5) and (A. 6) imply 



rf=Ki 



tjR, 



Near the border we can write 

oo 



(A.5) 

(A.6) 
(A.7) 

(A.8) 
(A.9) 



where the functions ,fn{tf) for n > 2 contribute to neither the metric, nor 
it's derivative at rj — Ri. Thus, we may freely set them to zero, fz{tf) = 
0, f/i{tf) =0, ... . Continuity of the metric at the border determine function 
fi{tf), while continuity of the derivatives of the metric at the border determine 
function /2(t/)- A similar procedure applies to tg. We expand it in a series near 
the bordering surface 



is=5o(V) + E(''/-^i)"5n(i/) 



The partial derivatives are given by 

drs 
drf 

dU 
dtf 

dts^ 

9rf r,=R, 



r!=Rl 



r!=Rl 



rt=Ri 



-^tyRi 
.fi(tf) 

9o{tf) 
9i{tf) 



(A.IO) 

(A.ll) 
(A.12) 
(A.13) 
(A.14) 
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Inserting these into Eqs.(A.2), (A. 3) and (A.4| we get 



500 



.901 



1 - 



1 4-2 



= 1 



rf = Rl 



910 = 



511 



9 \ „2 



rf=Ri 



rf=Rl 



Eqs.(A.8), (A.15) and (A.17) yield 



and 



50 (i/) 



1 



tlRi~r„ 



Ri{t]Ri-rg + -Rf) 



hits) = —^t}R,itjR, - rg)+gf{tf){tjR, - r. f 
t^Ri 



which imphes according to Eq.(A.16l 

9i{tf) 

hence 



CIV 



i{t)Ri-rg) 



1 



tjRl(tfRl-,g 



(A.15) 



(A.16) 
(A.17) 



(A.18) 



(A.19) 



(A.20) 



(A.21) 



Matching conditions follow from continuity of the derivatives of the metric, 
namely 



00 



dgi 

drf 



dgoi 
drf 



rf=Rl 



rf = Ri 

^ d9io 
drf 







d9_ 
drf 



11 



= 



(A.22) 



(A.23) 



(A.24) 



rf = Rl 
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Eq.(A.22| can be written in the form 

2 



dgoo 
drf 



r;=Ri 



rg dr,, ( dt^ 



r1 drf \dtf 



+ 1 



/ dr 



(r, - rg)2 drf \dtf 



rs J dtf dtfdrf 
1 ^dr, 



I /i 5o' + ( 1 



2 90 91 



dtf dtfdrf 
rn . 4 



rf = Rl 



(rs - rgY 



h-^tf-^Rl 



1 



2-tf-'R^h 



= 



(A.25) 



where 



h{tf) = 
9i{tf) = 



2tjRi-rg + lRl 



it}^Ri{t}Ri 
&Rl{t}Ri~r, 



9{tjR,-rg)^ 



(A.26) 



(A.27) 



according to Eqs.(A.20l, (A.21). Hence the Schwarzschild radius is given by 

4 



from Eq.(A.25|. The partial derivatives are given by 

/i(i/)=i| 



drs 
drf 



dt^ 
dtf 

dU 
drf 



ri=Rl 



9o{tf) 

9l{tf) 



9t] 



9tl 



4i?2 



6i?it/ 



9t] 



4i?2 



(A.28) 

(A.29) 

(A.30) 
(A.31) 



Integrating Eq.(A.30) we get 



Ri 



artanh t 



f ' 



iRi 
rn 



(A.32) 



therefore 



ts = tf + 3—tf — 3 
i?i 



Ri 



artanh t 



f ' 



J2irf-Rir9nitf) (A.33) 
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where functions gn{tf) for n > 2 - like in case of the coordinate ~ do not 
contribute at r/ = Ri. Function gi{tf) is aheady known. Functions /2(i/) and 
g2{tf) can be derived from Eqs.(A.23), ( |A.24 ), because 



dgoi 



drt 



rt=Ri 



dg 
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drf 



rj=R 



fi go gi 



{rs - rgY 



fi /o /i 



1 - 



1 ^ „ ) (-91 91 + .90 2 32) 
(/i /i + /o 2 /2) 

(A.34) 



dgi 



drf 



rf=Rl 



/i 5? + 1 



2 51 2 52 



1 



7 ^ ^ fl fl — FT 2/12/2 



1 



, 



hence 



and 



52 (t/) 



3tf (9tj - 8i?f 
(9t|-4i??)' 



(A.35) 
(A.36) 

(A.37) 



Finally we get the coordinate transformation between the flat FRW and the 
Schwarzschild metric as follows 



^9 A 



tf + i^n -3 



'a 

Ri 



artanh t 



+ (ry-i?i) 



9i| - ARj 



9t] - ml 



and 



together with 



=^|ry + (ry-i?l)■ 



(A.38) 

(A.39) 

(A.40) 
(A.41) 



Note that throughout we used units ensuring |i7o = 1- Return to general units 
is possible by the substitutions tf — > tf/to and i?i — > aoRi- 
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B Matching the closed FRW metric to the Schwarzschild 
metric 

The metric at the bordering surface is given by 



500 
9oi 
911 

922 
933 
902 



1 



1 - 



9lQ 



dtr 

ra = R2 

1 Sr., Sr., 



1 - 



J dtc drc 1 ~ ^ (^tc drc 

2 1 / n \ 2' 



r„=_R2 



-r1 Sir? d\ 



1 



1 



= 



920 — 903 ~ 930 = 



furthermore the derivatives at the surface are 
^500 



dgoi 
drc 

d9ii 



= 



ra = R2 



rc=R2 



dgio 
drc 



ra = R2 



drc 



-2a{tc) 



^ 



i?2 



rc=R2 



The coordinates can be written in the foUowing form 

rs = foix) + fi{x){rc-sma) + f2{x){rc~smaf 
ts = go(x)+9i{x){rc-sma)+g2{x){rc-siiiaf 

fo ~ P sin^ X sin « 



where 



according to Eq.(B.4). From Eqs. (3.6| and (B.l) we get 

— = 2/3- f-Wsm acot^x + l--r 

dx l-fj fo 



and from Eqs.(B.2) and (B.3) we have 



/i = /g^Jsin^acot^x+l-J 

Jo 



9i 



cos a 
tan a sin x cos x 



1 



(B.l 
(B.2 

(B.3 

(B.4 
(B.5 
(B.6 

(B.7 
(B.8 
(B.9 



(B.IO 
(B.ll 



(B.12 
(B.13 

(B.14 
(B.15 



fo 
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Hence we can express the Schwarzschild radius from Eq.(B.7) as 



Tg — (3 sin a 



Integrating Eq.(B.13) we get 

(3 cos a (l + 2 sin^ a) x — cos a sin x cos x 



90 



-2 sin a artanh 



tanx 
tana 



and 



91 ^ 13 



i sin^ X 

tan a sin'^ x cos x 
cos^ a — cos^ X 



From Eqs.(B.8l and (B.9) we have 

, 3 sin a 



/2 = /? 



4 cos^ a 



^Isin xcosx(cos a — cos x) 

g2 = p ^ 

2 cos'^ a (cos^ a — cos^ x) 



(B.16) 



(B.17) 

(B.18) 
(B.19) 

(B.20) 
(B.21) 



Hence the coordinate transformation at the bordering surface between the closed 
FRW and the Schwarzschild metric is the following 

ts = /3^cosq;(1 + 2sin^ a)x — cosasinxcosx 

o / tan x\ / ■ \ tan a sin^ x cos x 

— 2 sm a artanh ] + (i"c — sm a) ^ 

\ tan Of / """^^ ■ 



cos^ a — cos^ X 



1 2 sin X cos X (cos a — cos^ x) 
+ - (re — sin a) ^ — 

2 cos^ a (cos^ a — cos^ x) 



(B.22) 



n I -2 3 , . , 2 sin a 

rs^p[rcSin x+7(fc-smaj ^ 

' 4 cos"^ a 



together with 



(B.23) 



(B.24) 
(B.25) 



References 

[1] A. Einstein, and E. G. Straus, Reviews of Modern Physics, 17, 120 (1945). 



36 



[2] E. V. Linder, Astron. Astrophys. 326, 23 (1997). 

[3] V. Marra, E. W. Kolb, S. Matarrese and A. Riotto, arXiv:astro- 
ph/0708.3622v3 (2007). 

[4] V. Marra, E. W. Kolb and S. Matarrese, Phys.Rev. D77, 023003 (2008). 

[5] R. Kantowski, Phys.Rev. D68, 123516 (2003). 

[6] R. Kantowski, |arXiv:astro-ph/03084r9| ^l (2003). 

[7] T. Bisvas and A. Notari, [arXiv:astro- ph /0702555| >^l (2007). 

[8] K. Bolejko, Gen.Rel.Grav.41, 1737 (2009). 

[9] J. R. Mureika and C. C. Dyer, Gen.Rel.Grav.36, 151 (2004). 

[10] O. Gron and P. D. Rippis , Gen.Rcl.Grav.35, 2189 (2003). 

[11] K. Enqvist, Gen.Rel.Grav.40, 451 (2008). 

[12] M. Serpico, R. D'Abrusco, G. Longo and C. Stornaiolo, Gen.Rel.Grav.39, 
1551 (2007). 

[13] Sh. Khosravi, E. Kourkchi, R. Mansouri and Y. Akrami, Gen.Rel.Grav.40, 
1047 (2008). 

[14] C. Sormani, Geom.funct.anal.l4, 853 (2004). 



37 



